Frustrated total internal reflection and the illusion of superluminal propagation 



Vera L. BrudnjO 
Departamento de Fisica, Facultad de C'iencias Exactas y Naturales, 
Universidad de Buenos Aires, Cd. Universitaria, 
Pabellon 1, C14-28 EHA Buenos Aires, Argentina 

W. Luis Mochai|t| 

Instituto de Ciencias Fisicas, Universidad Nacional Autonoma de Mexico, 
Apartado Postal 48-3, 62251 Cuernavaca, Morelos, Mexico 
(Dated: July 15, 2008) 

We analyze the propagation of a pulse across a vacuum gap separating opposite flat parallel faces 
of two transparent dielectrics by means of an explicitly causal and retarded propagator constructed 
directly from the free-space wave equation. Nevertheless, our approach yields apparently superlu- 
minal propagation for the case of frustrated total internal reflection (FTIR) , that is, a transmitted 
wave packet appears on the far side of the gap at the same time that the corresponding incident 
packet crosses the front one. Thus, in this example superluminality is just an illusion, being con- 
sistent with both casuality and classical electrodynamics. We study the origin of the apparent 
superluminality in this case, which is inherent to light pulse propagation in free space and does not 
depend on the particulars of light-matter interaction, and find that it is due to propagation from 
the lateral wings of the incident pulse to the central part of the transmitted pulse. Thus, notwith- 
standing their similarities, FTIR is not equivalent to ID tunneling. We propose experiments to test 
our explanation of superluminality using opaque screens to block part of the wavefront, although 
we demonstrate that the propagation of smooth finite pulses constrained to be made up completely 
of evanescent Fourier components is indistinguishable from truly superluminal propagation, i.e., it 
may be completely accounted for using an explicitely superluminal and acausal propagator as well 
as the causal subluminal one. 

PACS numbers: 42.25.Bs, 41.20.Jb, 03.65.Xp, 42.25. Gy 



I. INTRODUCTION 

The physics of light propagation is a topic of active research, due to its relevance to several technological applica- 
tions and basic research. Recent research on new materials has shown that it is possible to exercise an extraordinary 
control on the propagation of light pulses, which in turn has generated new interest in both practical and fundamental 
questions on light propagation. Moreover, despite the fact that Maxwell completed the formulation of the classical 
theory of electromagnetism in 1864 and Einstein's special theory of relativity was presented in 1905, some hot con- 
troversies continue to arise on the subject of superluminal propagation [l|, |3, 0, M, H, S B S and its possible 
implications for both classical and quantum information theory pHIT^. Recent reviews on the subject can be found 
inRefs. 

One consequence of the special theory of relativity is that no signal can cause an effect outside the light cone of its 
source. Violation of this principle of relativistic causality leads to paradoxes such as that of an effect preceding its 
cause (l5l. [l6j. When dealing with light propagation in a material characterized by a given dispersion relation uj{k) 
between the frequency lu and the wave vector fc, several velocities may be defined, such as the phase velocity = uj/k 
and the group velocity Vg = V^rw. It is recognized that under certain conditions fl7'| both velocities can exceed the 
speed of light in vacuum c. This does not contradict the postulates of the special theory of relativity, for it has been 
recognized since the works of Sommerfeld [l^ and Brillouin [ll| that in order not to violate the principle of causality 
it is the information velocity that must not exceed c. The question that arises then is what is an appropriate and 
operative definition of information velocity. There have been several discussions and proposals on this subject, but 
the question is not yet settled [13, [U, [H 123, HI HI ■ 

The vast majority of the published work concerning superluminal pulse propagation deals with light propagation in 
material media, and the usual analysis attempts to explain how the interaction of electromagnetic radiation with the 
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medium affects the propagation of light pulses in such a way that they appear to travel superluminally. In a pioneering 
paper [5^ , Icsevgi and Lamb performed a theoretical investigation of the propagation of intense laser pulses through 
a laser amplifier. Apparent super luminal light propagation has been reported in gain-assisted systems [13, [IE US] 
as well as in birefringent crystals [13, [U, [13] , composite media and photonic crystals [H, [13] and dispersive media 

MM- 

There have also been claims of evidence of superluminal propagation in free space [371 |38| and during optical 
tunneling in frustrated total internal reflection (FTIR) configurations ^W, '4l', '42']. Optical tunneling has been 
studied by several authors, both theoretically and experimentally 44, 45, 46, 47]. Some of the observed results 
are still subject of debatable interpretation [13] and do not close the subject of whether there are possibilities for 
superluminal transmission of information in such systems. We have therefore chosen to address this subject in way 
that leads to straightforward interpretation of the results while resorting only to classical electromagnetic theory. 
We claim that although some results may appear to indicate superluminal propagation, there is no real superluminal 
transfer of information. 

In this paper we analyze mathematically the propagation of a pulse across a vacuum gap separating opposite flat 
parallel faces of two transparent dielectrics by means of an explicitly causal and retarded propagator constructed 
directly from the free-space wave equation. Our results yield indeed an apparent superluminal propagation corre- 
sponding to the conditions of FTIR, but they show explicitly that it is consistent with both casuality and with classical 
electrodynamics. Our example shows supcrluminality effects inherent to light pulse propagation in free space which 
therefore does not depend on the particulars of light-matter interaction. The illusion of supcrluminality consists 
of transmitted pulses arriving to the far side of the gap in synchrony with the crossing of the front surface by the 
incident pulse. We explain this illusion of superluminal behavior in terms of a causal, subluminal propagation, taking 
into account the spatial extent of the incident pulse along its transverse as well as its longitudinal directions and 
we propose experiments to demonstrate the retarded nature of propagation in FTIR. Nevertheless, we find that for 
constrained pulses fully made up of evanescent Fourier components, subluminal and superluminal propagation in 
FTIR experiments are indistinguishable. 

This paper is organized as follows. We first introduce a propagator that will allow the analyzis of electromagnetic 
pulse propagation across a vacuum gap (Section |TT|. This propagator is both casual and retarded and complies with 
the classical electromagnetic theory. We then study the propagation across the gap of pulses that arrive as plane 
waves with well defined angles of incidence. In Section IIIII we study the case of subcritical angles, yielding non- 
evanescent transmitted waves. In Section HVl we study the case of hypercritical angles, yielding evanescent waves. We 
conclude that the propagation of a light pulse in a FTIR configuration may appear superluminal and acausal but 
that it is actually subluminal and that propagation has to account necessarily for the lateral wings of the incident 
pulse. In Section |V] we suggest experiments that might demonstrate the actual causal and subluminal nature of 
the apparent superluminal behavior by using sharp opaque screens that block parts of the incident wave so that its 
extent becomes finite along both its propagation and its transverse directions. Nevertheless, as the borders of these 
screens produce propagating diffracted waves, in Section |Vl] we eliminate them and we study incident pulses that 
are finite along several spatial directions but that have a smooth profile. We obtain that if they are comprised of 
evanescent Fourier components only, they appear to propagate superluminally through the vacuum gap, even though 
their behavior is determined by our causal retarded propagator. In Section IVIII we construct an alternative acausal, 
superluminal propagator, and prove that it is exactly equivalent to the causal and retarded propagator when applied to 
fully evanescent finite pulses. Thus, for such constrained pulses, it is impossible to distinguish causal subluminal from 
acausal and superluminal propagation; the illusion of supcrluminality appears to be not only a matter of interpretation 
of the result of the propagation, but may be also present in the description of the propagation process itself. We present 
our conclusions in Section IVIIII 



II. PROPAGATOR 



Consider two transparent dielectrics occupying the regions z < and z > d. In this section we obtain the causal and 
retarded propagator that describes the motion of a pulse across a vacuum gap < z < d spanning from the interface 
at z = to that aX z — d (Fig. [T]). As we want to consider explicitly the angle of incidence onto the interface z = 0, 
we cannot treat the problem beforehand as if it were ID. For simplicity, we will assume full translational symmetry 
along the y direction, so that our problem becomes 2D. Thus, we start with the scalar wave equation 

(V^ - l^)Go(t, X, z- 1', x', z') = 5{t - t')S{x - x')6iz - z') (1) 



' d ' 

FIG. 1: Pulse with arbitrary profile incident at an angle 9i < 9c onto the surface z — of a, dielectric with index of refraction 
n = 2. The pulse is transmitted into a vacuum gap and into a second dielectric at z — d. The speed of propagation within each 
media, c/n and c, are indicated, as well as the speed of propagation of the wavefronts along the interfaces 



with a unit singular point source fired at time t' at position (x',z'), which is solved by the explicitly causal and 
retarded free space Green's function j48|] 

^ e Ut - 1') - ^ [x - x'Y + {z - z'y\ 

Go(t,x, z't' ,x' , z') = , , (2) 

^ ^ 2^ ^c^t - t')^ - {x - x')^ - {z - z')^ 

where Q{. . .) is the Heaviside unit step function. Using image theory we can construct a Green's function that obeys 
Dirichlet boundary conditions at the surface z — 0, 

G{t, X, z] t' , x' , z') = Go{t, X, z; t' , x' , z') — Go(t, x, z; t' , x\ —z'). (3) 

In the half-space z > Green's theorem yields the solution 

= P(,,,i',.'.0-)*(*'.y.O") (4) 

of the homogeneous scalar wave equation that is outgoing as z ^ oo and is null in the remote past, where x', 0^) 
denotes its previous values on the boundary z = 0^ , and 

P{t, X, z- 1', x\ 0+) = ^ G{t, X, z- 1\ x', z')L,_o+ (5) 
az' 

is the causal retarded propagator of the problem. We may identify the field (f> with the component Ey of the electric 
field E in the case of a TE or s polarized incoming wave, and with the component By of the magnetic field B in the 
case of TM or p polarization. The propagator ([5]) does not account for the presence of the two dielectrics bounding the 
air gap. Thus, the field ([4]) has no information about the multiple reflections at the boundary of the gap. In principle, 
these can be incorporated by reflecting the fleld at the interfaces z = 0, d using the appropriate Fresnel coefficient and 
propagating it back and forth across the gap with the propagator ([5]) for the z — Q surface and a similar one for the 
z = d surface. The total field would then be the sum of all the multiply reflected fields and would have information 
about the electromagnetic properties of the reflecting surfaces. In this paper we will restrict ourselves to an analysis 
of the first crossing of the air gap 0^ — s- , and thus our results will be unrelated to the nature of the bounding 
media. 
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Substituting Eq. ^ into ([3]) and ([5]) we obtain 



^ g Q (c{t - t') - ^{x~x'Y + zA 

P{t,x,z-t',x',Q+)=--^ \ A. (6) 

TT az ^c2(t-i')2-(x-a;')2-z2 



The field (p can then be written in terms of an ancillary function 



where 



d 

(l){t,x,z) = -—il){t,x,z), (7) 



^{t, x,z) = - f dx' I dt' Ht',x',0+) 

TT J J ^c^{t-t'y-{x-x'y-z^ 



plays the role of a potential and the integration is performed within the region c{t — t') > ■^/(a; — x')^ + z^. Clearly, the 
procedure above yields a causal (sub)luminal propagation from the z = 0+ plane to any point in the z > vacuum. 



III. NON-EVANESCENT WAVE TRANSMISSION 



We consider now that an arbitrarily shaped pulse impinges at a well defined angle 9i on the inside surface z = 0^ of 
a homogeneous non-dispersive dielectric with index of refraction n (Fig. [1]). The incident pulse is therefore described 
by an arbitrary function (j)i{t, x, z) = fi[t— in/ c)ni • p] of a single variable t — {n/cjfii ■ p, where fii = (sin 9i, cos Oi) is 
a unit vector pointing along the angle of incidence 9i and p = (x, z). Notice that in this case the incident wavefronts 
have an infinite extension in the direction normal to hi. At z = 0"'' and after being transmitted into vacuum, the 
outgoing field can therefore be written as 

(l)it,x,0+) = ftit-x/v\i) (9) 

where ft is related to the arbitrary function fi and the Fresnel amplitude for transmission from the dielectric into 
vacuum. As we are concerned only with propagation across the vacuum gap, we will take ft as given and we will 
disregard its relation with /i, which would involve the dielectric properties of the incident medium. 

The intersection of the pulse with the interface z = is therefore seen to travel along x with velocity W|| = c/{n sin 0i) 
(not to be confused with the parallel component of the incident velocity hic/n). For incidence angles smaller than 
the critical angle 9 < 6c = sin^^(l/n), v\\ > c and we have normal transmission, while for 9 > 9c, v\\ < c and total 
internal reflection ensues. 

In order to set up a reference with which to compare the evanescent case, in this section we employ our propagator 
to study a non-evanescent plane pulse. Thus, we consider here the case 9 < 9c and we substitute Eq. ^ into Eq. ([8]) 
to obtain 

iPit, x.,z)^^ j dt" ft{t") j dx' {c\t - t" - x'/v\\f ~{x- x'f - z'^Y^''' (10) 

after changing integration variables from t' to t" = t' — a;'/^||- The integration (fTO|) has to be performed over the 
region c{t — t") — px' > \J (x — x'Y + z^, where ^ = c/t;||. Thus, t" has an upper bound 

tm =t - {px - VZ)/C, (11) 



where v = ^\ — c2/w|, and for each value of t" <tm, x' is bounded by the limits 

X'^ = -(l/l^') (x - pc(t - t") ± ^[px - C{t - t")Y - I^2z2) . (12) 

A simple interpretation of Eqs. pT|) and p2)) can be obtained with the help of Fig. [2l Consider an event £ = {ct, x, z) 
defined by the observation of the field at a given position (x, z) with z > and at a given time t. Causality requires 
that only events within the past light-cone of £ are able to influence it. Notice that a given value of t" denotes a 
point that moves along x keeping a fixed position with respect to the intersection of the incident pulse with the z = 
interface. In Fig. [2] we show the world lines {ct\x',0) of a few such points. Since each of them moves with speed 
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FIG. 2: Contributions to the non-evanescent wave observed at an event £ = {t,x,z). Some world lines {ct',x' = — i")iO) 
of source points labeled by fixed values of t" and thus moving on the z = plane along with the incident pulse are shown. 
The positions x'j- denote the intersections of the past light cone of £ with one of the world lines plotted. We consider the case 
V|l > c. The coordinate and time axes are indicated. 



W|| > c, most of its world line lies outside the past light cone of £. Only if t" < t„i can it actually cross the light cone, 
entering and leaving at positions x'_ and x'j^ respectively. 
After the change of variables from x' to 



Eq. (fTU]) simplifies to 



_ x' - [x - fic{t - t")] 
V = — , (Id) 

^[nx-cit~t")]^ -i^^z^ 

^{t, x,z) = ^ /*" dt" Mt") f -^M=. (14) 

J-oo J-1 ^/l - if 

The integration over 77 is immediate, so that substituting Eq. ([14] in Eq.((7]) we obtain finally 

(j){t,x,z)^ ft{t-hf p/c). (15) 

where fit — (/i, v). Thus, as illustrated in Fig. [U the transmitted wave is a pulse with the same profile as the incident 
field and propagating with speed c at the well defined angle 9t = sin^^ /i = cos^^ z^, in accordance with Snell's law as 
could have been expected. 



IV. EVANESCENT WAVE TRANSMISSION 



We consider now the case 9i > 9c, for which U|| < c and the transmitted wave becomes evanescent. In this case we 
substitute Eq. ([9|) into Eq. ([8]) to obtain 

ij{t,x,z)^^ J dt"ft{t^x/v\i+t") J dx" {c'it" + x" /vi\f - ix"f - z^y^^\ (16) 

after introducing the variables x" = x' ~ x and t" = t' — t — x" /v\\. As shown in Fig. [3l for any observation event 
£ and any value of t", there is exactly one intersection a:" between the past light cone of £ and the world line 
(ct'j x" — W|| {t' — t" — t), 0), where now x'!_ — —^(3['yct" + ^/z^ + (7/3ci")^] and we introduced the definitions f} = v^j/c 
and 7 = 1/y^l — wj^/c^- We have assumed that U|| > 0. As the world line leaves the past Hght cone at x'L, in Eq. ([16]) 

the integration over t" is unconstrained and that over x" extends from — c» to a:" . Another change of variable, from 
x" to 

1 x" + -f^l3ct" 



7/3 + (7/3ct")2 ' 

yields 



(17) 



^(i, x, z) = ^ r dt" Mt ~ xlv\\ + t") [ ' -j^. (18) 
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FIG. 3: Contributions to the evanescent wave observed at an event £ — {ct,x,z) as in Fig. (2] but for iiy < c. The position x'_!_ 
denotes the intersections of the past light cone of £ with one of the world lines x" = — t" — t) corresponding to t" . 



The integrals over x" and over rj in Eqs. I|16p and (jlSp respectively yield an infinite value for ip. This divergence 
is not unlike that commonly found for the electromagnetic potentials produced by infinitely extended sources. For 
example, when calculating the electric field produced by a uniformly charged plane one cannot simply obtain the 
corresponding potential by integrating the Coulomb kernel over the whole surface. However, in that case the electric 
field may be obtained either by deriving the Coulomb kernel first and integrating afterwards or else, by truncating the 
integrations at a finite distance, deriving the resulting potential to obtain the field and afterwards taking the limit of 
an infinite surface. Here we follow the later procedure. Thus, we set a finite lower integration limit x'J^ in Eq. (jl6p . 
corresponding to a lowest point rjL in Eq. (|18p . and we take the limit x'[ — oo, rjL ^ ^oo after obtaining the field 



As tp depends on z only through 77^, substituting Eq. (fT5|) in ([7]) we obtain 



cp{t^x,z)^:^ 1°" dt"Mt-x/v^\+t")^^^. (19) 



TT 



In the limit x'[ —t —00 we evaluate 



- 7 r^__ -,2 , (^,R^i.ii\2' (2*^) 



and we obtain finally 



Taking the absolute value of v\\ in the numerator of Eq. I|2ip allows its use also for v\\ < 0. 
To grasp the meaning of Eq. (|2ip we evaluate it for an infinitely sharp pulse 

/t(T) = /o<5(r). (22) 

Substitution into Eq. ([2T]) yields 



nt, x,z)^ -2 . (23) 

TT + 7^(2; — W|| t) 

Surprisingly, at any time t the pulse transmitted at a distance z from the interface is given by a Lorentzian of width 
2/7 centered in front of the actual position a; = of the incident pulse on the z = surface (Fig. 3]). Thus, the 
propagation seems to be instantaneous in the direction normal to the surface, and actually, part of the pulse seems 
to travel backwards in time [39l. [47|. as at a position (x,z) it becomes appreciable at times t < s/wy, that is, before 
the incoming pulse reaches the corresponding position (x,0). However, our deduction of Eq. shows that it is 

completely consistent with a causal and retarded propagation, and that the field at a (x, z) at time t is not produced 
instantaneously by the incoming field at (x, 0), but arises from previously excited positions {x' , 0) with x' < x + x'!^. 



FIG. 4: Infinitely sharp wavefront (heavy solid line) incident at an angle 9 = 34.5° > 6c upon the surface of a dielectric with 
index of refraction n = 2. The pulse widens and diminishes as it is transmitted across an air gap and into a second dielectric. 
The thin lines indicate the nominal pulse width and the dashed line its center. 

It is interesting to note that, according to Eq. ([23]) . the height of the transmitted pulse is inversely proportional to 
the distance z from the surface, instead of decaying exponentially as usually found for evanescent waves. However, 
Eq. (j23p describes the propagation into vacuum of a single infinitely sharp incident wavefront. In a wavetrain made 
up of a succession of incident pulses, the regions excited by neighboring pulses overlap each other, as the width of each 
transmitted pulse increases in proportion to z, and therefore their corresponding fields interfere. This interference is 
at the origin of the exponential decay of periodic waves, as can be verified by choosing 

/t(r)=Ae-^"^ (24) 

and substituting into Eq. (pi]) . which yields 

A simple contour integration closing the integration path with an infinite semicircle on the lower half complex t" 
plane yields the familiar result 

0(i, X, z) = Ae'('3^-"*)-«^ (26) 

where Q = uj/v\\ = nsinOiCo/c is the parallel component of the wave vector and k = a;/(7|w|||) = \JQ^ ~ uP' j (? — l/l 
the inverse of the decay length I. As increases, the distance along x between successive maxima and minima 
decreases, yielding larger interference effects and a shorter decay length. 

V. SCREENS 

The results of the previous section suggest an experiment that could confirm that propagation of evanescent waves 
in the FTIR geometry is not superluminal nor acausal. The experiment could be performed simply by partially 
covering the surface of the first interface with a couple of opaque screens as shown in Fig. [5l If transmission were 
indeed superluminal, we would expect a non-null transmitted field (f){t,Xa,d) across the gap in front of the edge Xa 
of the first screen Sa as soon as the leading wavefront J-^ of the incident pulse reaches Xa- Similarly, we would expect 
that the field (j){t, Xf,, d) would be modified as soon as reaches the edge Xf, of the second screen St,- 

To calculate the field corresponding to Fig. [5] we go back to Eq. ((TB)l . The screens confine the integration region 
to the interval Xa < x' = x -\- x" < Xb- This inequality has to be obeyed together with the previous constriction 
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FIG. 5: Pulse as in Fig. [T]but impinging on the surface of the first dielectric at an angle 9i > 9c. The surface is covered by 
semi-infinite opaque screens Sa and Sb with edges at Xa and Xb- The leading wavefront of the incoming pulse is about to reach 
Xa- T denotes the foremost wavefront. 



x' < X ^ x"_. These conditions can only be satisfied by those wavefronts which have reached Xa and left behind the 
first screen before t — Ta, where Ti = (l/c)-\/ {x ~ Xi)^ + z"^, i ~ a,b, is the minimum time required to reach {x,z) 
from {xi, 0) moving at speed c. Thus, only those points on the wavefront labeled by t" > T" can contribute to (|16|) . 
where 



rpll 



— \/ {xi — x')^ + z^, i — a,b. 



(27) 



For those wavefronts which have left Sa by time Ta but have not reached Sb at time Tj,, namely, those with T^' < t" < 
T", the integral over x" in Eq. has to be performed from x'a = x + Xa up to . Finally, for those wavefronts 
which have already been blocked by Sb by time Tf,, namely, those with t" < , the upper limit of integration has to 
be replaced by x'^ = x + x^- Therefore, 



^{t, x,z)=^ r° dt" Mt - x/v^^ + 1") 

IT 7t/' 



dri 



(28) 



where we used the change of variables (fT7|) and substituted a;" — > a; + Xi in it to define the limits rji . Notice that ij) 
depends on z only through the integration limits rji, so that substituting Eq. (|28[) in ([7]) we obtain 



'ip{t^ X, z) 



7/3c 



where 



is given by Eq. ^and Q = C{x'-,t") with 



dt" ft{t - x/v^\ + t")Kae(T: - 1") - 6e(r," - 1'% 

x" + j^Pct" 



^/{x"y + 2-i^(3cx"t" - 72/32(^2 _ c2(t")2) 
Substituting Eqs. and in ([29]) we finally obtain 



(j){t,x,z) ^ - I dt" ft{t - x/vH +t") 



7Vj|x 



z2 + (7U||i")' 



r(l-C(n), 



(29) 



(30) 



(31) 



(32) 
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FIG. 6: Field isolines of (j>(t,x,d) produced by an unit delta function input pulse that propagates along the z = Q surface with 
speed = and is blocked at a; < a;^ = —3d and a: > it = 3d by opaque screens (shaded regions). The world line a; = nyt 

of the incident pulse is indicated by a dashed line. The times ta and tb when the incident pulse cross the edges of each screen 
are indicated by the horizontal dot-dashed lines. The singularities of the transmitted field are indicated by the thick solid 
hyperbolas. 



where 

c{t") = 1 + [C6e(r," - 1") - qmt: - <")]■ (33) 

Notice that the field in the presence of screens (Eq. (P^ ) is given by an expression similar to that corresponding 
to the field in the absence of screens (Eq. ([21])) but with a correction term C due to the diffraction by the screen. 

As in the previous section, we consider again the case of a sharp incident pulse Eq. (|22p . Substituting in (|32p we 
obtain 

{ if ct<Xa/l3+y ^{x-Xaf + Z^ , 

(pit, X, Z) = ^2 _^ 21^ _ y 1)2 ^ \ (i^a - X, X / - t) ~ (i^b - X, x/v\\ - t) if ct > Xb/(3 + {x - Xb)'^ + Z^, (34) 

z 7 [x t>|| ) ^ ((xa — x,x/v\\ — t) otherwise. 

We remark that the field is zero until the time Xa/v\\ when the incident pulse shows up from behind the screen Sa, 

and this information has had enough time (l/c)-y/ {xa — x)^ + z^ to propagate from the screen's edge {xa,0) to the 
observation point {x,z). Similarly, the information that the pulse has hidden behind screen Sb does not reach the 
observation point until the time Xb/wy + (l/c)\/(x — Xb)^ + z^. The field has singularities due to the passage of the 
incident pulse through the screen edges, that propagate at speed c from the events {xa/P,Xa,0) and {xb/P,Xb,0)- 
Notice that (^(xi — cc, 0) ^ 1 as a: — *■ oo. Thus, if we follow the incident pulse, i.e., we take x « v\\t, then cj) 
asymptotically after the pulse hides behind Sb- Furthermore, if the screens are very far apart we recover the field (j23|) 
between the screens. 

The features above are illustrated in Fig. [5] which show the transmitted field at the plane z = d. Notice the delay 
d/c after the incident pulse crosses Xa at ta = Xa/v\\ before a non-null field first appears across the gap at (xa, d), and 
a similar delay after the incident pulse crosses Xb at tf, — Xf,/u|| before the field starts to be extinguished at (xb,d). 
Furthermore, notice that for some time the field penetrates a small distance « d beyond Xf, as if there were no screen. 
The field is singular at the hyperbolas with vertices at x — Xa, t = ta + d/c and at x = Xb, t = tb + d/c given by 
the intersection of the {ct,x,d) hyperplane and the future light cone of the events {cta,Xa,0), {ctb,Xb,0)- Thus, we 
have shown observable consequences of the fact that evanescent waves in FTIR do not propagate superluminally nor 
acausally in the direction normal to the dielectric-vacuum interfaces, but with retardation and obliquely. A graphical 
approach to the results of this section and an animation illustrating them may be found in Refs. |40 and Ref. fisi 
respectively. 
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VI. A SMOOTH TRANSVERSE PROFILE 



Perfectly opaque screens such as those considered in the previous section introduce sharp discontinuities in the pulse 
at z = 0. The truncated pulse no longer has a well defined propagation direction 9i but may still be represented by a 
superposition of pulses with varying propagation directions. A sharp truncation leads to the presence of subcritical 
incident angles di < 9c, and therefore to the presence of both, evanescent and non-evanescent transmitted fields. It has 
been argued [41| that the retardation effects discussed in the previous section may be due only to the non-evanescent 
contributions, known to be subluminal. The comparatively slow subluminal contributions would be unable to affect 
the arrival of the superluminal signals if the later were actually present. However, any small non-evanescent wave 
would dominate the transmitted signal after a wide enough gap. Thus, it is interesting to study the propagation of 
pulses with a finite transverse extension but with a smooth lateral cutoff and built up completely from hypercritical 
9i > 9c evanescent contributions. 

To explore the propagation of the smoothly truncated pulses discussed above, we consider an incoming field given 
by a Fourier integral 

+ /"do; f dQ ,(Q^^^t) 



0(<,x,O+)= j — J ^/^ge^W--*), (35) 

where fi^g is the amplitude for each parallel component of the wave vector Q and frequency lo. We can change 
integration variable from Q to the parallel velocity v = lu/Q, 

<^(t,a:,0+)=.|g||/..e-(-^/''). (36) 

where we introduced the velocity dependent amplitude = {i^/v'^)fi^ i^/y. The incident field (j36p will give rise to 
evanescent waves exclusively as long as all non-null components f^jv have v < c. At this point we could integrate first 
Eq. with respect to uj, obtaining thus a superposition of plane pulses, each of which may be propagated across 
the air gap according to Eq. (j21[) . Alternatively, we may propagate each monochromatic component using Eq. (|26|) 
and afterward perform the integrations in Eq. (|36p . We follow the later approach and write 

0(t,x,^) = / ^ / (37) 



where j — Xj ^\ — j (? as in Sec. IIVI For simplicity we assume that we may factor /^^i, — fujfv into frequency and 
velocity dependent amplitudes, fi^ and /„ respectively; controls the time duration of the pulse, or equivalently, its 
longitudinal extent, while /„ controls its transverse extent. We further assume a narrow Gaussian velocity distribution 
of width Av around a nominal velocity vq < c, 

= (38) 
Av 

where u = v — vq. We assume Av is small enough that the non-evanescent contributions to the field may be neglected 
and the exponent in Eq. (I37p may be linearized in u. Thus, the transmitted field becomes 

(t){t,X,z)K [ ^ f^f.~'^[^/ilovo)+i{t-x/vo)] f '^^ ^-u^ l2i^'X?+^u\zho-ix\lvl ^ j^gg^ 



where 70 = 1/ ■\/l — Ug/c^. The integration over u is immediate and yields 

(t){t,X,z) = y — f^Q-'-'l^/i'(ova)+i{t-x/vo)]-[Avuix+iz/^„)f /2v^ _ ^^q-j 

In Fig. [7] we illustrate the results of applying Eq. (HO)) to a pulse with a Gaussian frequency distribution, 

f^^A^e'^^-^"^"/^^^", (41) 
Aoj 

of area A and width Aw centered at luq {ujo — 16c/d, Aco = 2c/c?, vq — 0.7c, Av = 0.15c). The pulse is seen to form 
on the z = surface at a; ~ —5d and at time t « 7d/c, it propagates along the x w v^t line for a while, peaks at a; = 
at time t = and disappears at a; « 5d, t « 7fi/c. Its maximum duration r and size L are ct k L k d for a fixed 
position and fixed observation time respectively, and it contains altogether about six nodes. 
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FIG. 7: Intensity of a Gaussian pulse localized in both space and time {x and t) incident on the front face of an air gap of 
width d (left) and after after crossing it (right). The nominal velocity along xis vo = 0.7c with width Av = 0.15c. The nominal 
frequency is ojq = 16c/d with width Atj = 2c/ d 



Surprisingly, after crossing the gap, the pulse looks essentially the same! It appears at the back z = d face of the 
air gap at roughly the same time and the same position as at the front z — face. It also peaks at the origin at 
t — and disappears from the back surface in concordance to the incident pulse on the front face. Thus, it truly 
appears to propagate instantaneously. The main difference between the incident and transmitted pulse is that the 
intensity of the later is suppressed by 12 orders of magnitude. Another interesting difference is that the number of 
visible nodes in the transmitted pulse has decreased to about 4. This is a consequence of the fact that in FTIR, 
the Fourier components with higher frequencies are more damped than those with lower frequencies. Finally, a more 
subtle difference is that the speed of propagation along the back face is slightly but noticeably larger than that on 
the front face. This is due to the fact that plane waves incident at angles closer to 9c have a larger penetration length 
than waves incident at larger angles. Thus, the angle of propagation of the transmitted pulse is smaller than that of 
the incident wave 

In Sec. II VI we have argued that the transmission of evanescent plane pulses across an air gap under FTIR conditions 
is fully consistent with a retarded and causal propagation along oblique directions. We have strengthened our argument 
by showing that there is a delay before a perturbation, such as blocking part of the incident wavefront, can produce 
an effect on the pulse transmitted across the gap. Furthermore, by truncating an incident plane pulse producing 
an abrupt transverse profile, we showed that the transmitted pulse is shifted along the surface in the direction of 
propagation. However, in Sec. IVII we showed through an example that if the pulse has a smooth transverse profile, 
such that all its Fourier components are evanescent, it is transmitted as if it were indeed super luminal. To understand 
this result, in Fig. [8] we show schematically a pulse smoothly truncated along its transverse direction, built from 
narrow plane components propagating along well defined directions 9 > 9c- The incident field peaks at the regions 
where the directional components add in phase, indicated in the figure by elliptical regions around the crossing point 
of the different incident wavefronts. As each of the components crosses the gap, it is widened according to Eq. (|23p . 
The peak of the transmitted pulse appears in the regions of largest overlap between the different widened transmitted 
components, i.e., at the regions where their centers coincide. The peak of the transmitted pulse at the back face of the 
air gap is seen to appear at the same time as the peak of the incident pulse reaches the front surface. Thus, it would 
seem as if the peak tunneled instantaneously in the direction normal to the gap. However, in the previous sections we 
have shown that each of the transmitted components of the transmitted wave originates causally from regions in the 
lateral wings of the incident pulse. Thus, the peak of the transmitted field does not actually originate from the peak 
of the incident field; it is formed by contributions from the lateral wings of the incident field, which reach the front 
face of the air gap first. The lateral wings of each component have enough time to cross the gap traveling at speed 
c and combine to form the relatively small transmitted peak right at the time when the larger incident peak reaches 
the front surface. Similarly, the different components of the field produced by the peak of the incident pulse get out 
of step as they cross the air gap and, therefore, do not contribute to the peak of the transmitted pulse, but rather, 
to its lateral wing. This is illustrated by the dashed arrows in Fig. [8l Thus, it seems that the physical propagation 
(retarded and causal) can not be distinguished from the nonphysical propagation (superluminal and non-causal) as 
long as we only consider smooth incident pulses which contain only propagation directions above the critical angle 
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FIG. 8: Evanescent transmission of a smoothly truncated pulse made up initially of a superposition of narrow plane wavefronts 
(heavy solid lines) with a distribution of angles above Oc- As it crosses the gap, each wavefront widens (shaded bands around 
dashed lines, as in Fig. |4|. The ellipses represent schematically the contour levels of the pulse and are centered at the region 
where the different contributions add coherently in phase, i.e., the regions where the centers of each component coincide. 
Snapshots are taken at three different times. The solid arrows illustrate the nominal propagation of the peak. The heavy 
arrow indicates the apparent instantaneous tunneling across the air gap. The dashed arrows illustrate the actual (sub)luminal 
propagation. 



MM- 

VII. EVANESCENT PROPAGATOR 



In Sec. IVII we found that a particular pulse made up of only evanescent components seemed to propagate instan- 
taneously across the air gap, in contrast to the abruptly truncated pulses considered in Sec. |V]for which retardation 
effects have observable consequences. This was explained graphically in Fig. [5]for incident pulses built up from narrow 
directional components, each of which is widened as it is transmitted across the gap. To show that this behavior is 
generic, we start from the Fourier decomposition of an arbitrary field 



(b{t,x,z) = J _eHW--')0^^(z), (42) 

where 

0^,q(0)= [ dt [ dxe-''-^^-'^'U{t,x,z). (43) 



The condition that (j)(t, x, z) is made up exclusively of evanescent waves is equivalent to stating that the integration 
region in Eq. ((35|) is given by \w\ < |Q|c, i.e., (/)i^,q = if |a;| > \Q\c. Using Eq. ([26]) we propagate each Fourier 
component from z = 0+ to z > as 

c^u..,q{z) = e-''^0^.Q(O+), (44) 



where k = \JQ'^ — cj^/c^. Thus, we can combine Eqs. (|42|) . ([43|) and (|44| to obtain 

0(t,a:,z)= / dx / dt'F'(t,a;,z;t',x',0+)<?!)(t',a;',0+), (45) 
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where 



2tt / 2tt 



(46) 



and the integration region is given by |w| < \Q\c. Comparing Eq. (|45l) with (|4]) we find that P' is a propagator that 
can be used in the same way as the propagator P defined in Eq. ([6|) to find the vahie of the field at z > given 
its values at z = 0^, provided the field is built up of evanescent components only, as in the examples of the two 
previous sections. We remark that our original causal, retarded and subluminal propagator P was able to propagate 
any arbitrary outgoing field. However, by adding constrains to the field, we gain freedom in our choice of propagator, 
as we can chose arbitrarily its effect on fields that do not obey the constrain. Thus, if we can find any function P" 
such that 

dx' J dt' [P"{t, X, z; t', X, 0+) - P{t, X, z; t' , x, 0+)]/(t' - x /v) = 0, (47) 

for an arbitrary flat pulse / moving along x with any velocity — c < w < c, then we could employ P" instead of P 
to propagate an arbitrary evanescent pulse. P' above is just one of the many possible choices of a propagator for 
evanescent pulses. 

To proceed, we make a change of variable u — + Qv to write Eq. as 



-c 27r 



(48) 



where 7 — l/-\/l — w^/c^, £^=x — x',T = t — t' and we perform the integration over 



dv 



{z/^f{£,^vrf 
[{zhY + (e - vrf 



(49) 



Notice that P'(i, x, z; t' , x' , 0^) is symmetric under the interchange ^ <-> — ^ and also under the interchange r — r, 
i.e., P'{t,x,z;t',x') = P'{t,x' , z;t' ,x) = P'{t',x,z;t,x') = P' {t' , x' , z;t,x). Thus, the evanescent propagator P' is 
superluminal and acausal. 

To finish the calculation of P' we make another change of integration variable v = c sin a to write 



27r2 



""^^ , (z cos Q!)2 - (C - CT sin a)2 

da COSQ!r- -T5 — : 

[(zcosa)^ + (C — cr sma)"^J^ 



and we perform the integration 



-tt/2 



z log 



(cV^ - z^ - zs)2 - c\^e 

[c^t"^ — Z^ + ZsY — C^T^^^ 



-4s 



where p — \/ S,"^ + z"^ is the spatial distance from the source to the observation point and 
space-time interval. This expression may be simplified to 



z\og{\z + s\/\z-s\)/s^ if s2 >0, 

-2zarctan(|s|/z)/|s|3 - 2/|s|2 if s2 < q, 



(50) 



(51) 



c T is the squared 



(52) 



The evanescent propagator is displayed in Fig. [S] The figure shows explicitly the temporal and spatial symmetry, 
and thus the superluminality and non-causality of P' . Notice that P' has singularities at the projected light-lines 
X — x' = ±c{t — t') which converge at the origin x = x' , t ^ t' . Thus, propagation is largest for instantaneous 
propagation in the direction normal to the air gap. 

We have found that the propagation of evanescent pulses can be described with either the exact propagator P of 
the problem, which is causal, retarded and subluminal, or with an evanescent propagator P' which is superluminal 
and non-causal. Both yield exactly the same transmitted pulse when the incident pulse contains only evanescent 
components. Thus, it seems to be impossible to distinguish superluminal from subluminal propagation in experiments 
performed with purely evanescent pulses [41]. 

An explanation for the curious conclusion found in this section can be obtained by going back to Eq. (|42p which 
we rewrite as 



4'it,x,z) = / ^(l)Q{t,z). 



(53) 
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FIG. 9: Propagator P'{t, x, d; t' , x' , 0"*") normalized to c/d^ as a function of the displacement x — x' and the delay t — t'. Distance 
is measured in units of d and time in units of d/c. For aid in visualization, the height of the propagator was truncated at 
P' = 0.1c/ d"". 

Notice that for each finite wave vector Q, the time dependent Fourier component 4>Q{t,z) has a strictly finite spectrum 
— |Q|c < oj < \Q\c. Thus, (/>Q(t, z) is an analytical function of t with no singularities and can be analytically continued 
to an arbitrary time t2 from its values in an arbitrarily small neighborhood of any other arbitrary time ti. Therefore, 
(j)Q{t,z) is perfectly predictable in principle. Using antropomorphic language, we may say that at time ti the system 
knows from the present values of (t>{t « ti,x, 0+) what its future values ({){t m t2 > ti + z/c, x, 0+) will be, and thus 
it can use that knowledge to build subhiminally a transmitted pulse (p{t w 1,2, x.z) that will mimic s=; 1,2- x,0^), 
giving the impression that super luminal transmission has taken place at time t2 ■ The validity of the argument above 
in the presence of thermal or quantum noise has to be investigated. 

VIII. CONCLUSIONS 

To study the propagation of a light pulse through a vacuum gap between two parallel dielectrics in a FTIR 
configuration we constructed a propagator derived directly from the wave equation resulting from Maxwell's equations. 
This propagator is retarded and complies with the relativistic causality principle inherent to classical electromagnetism. 
Therefore, it can only account for causal, (sub)luminal propagation of light pulses. However, when this propagator 
is used to study the propagation of wave packets through the gap, we find apparent superluminal behavior, that is, 
a wave packet might appear on the far side of the gap at the same time that the incident packet reaches the front 
one. Therefore this illusion of superluminality, present within classical electromagnetic theory even in vacuum, is fully 
consistent with relativistic causality. We showed explicitly that propagation in FTIR actually takes place subluminally 
between the lateral wings of the incident pulse and the central peak of the transmitted pulse, and we proposed simple 
experiments that could verify this statement. Thus, although FTIR has many similitudes to ID tunneling, its correct 
physical interpretation requires a 2D or 3D analysis. On the other hand, we constructed an explicitly superluminal 
and acausal propagator that yields identical results as the retarded causal one when applied to smooth pulses made up 
of evanescent contributions only. Thus, there is a class of pulses for which superluminal and subluminal propagation 
would be indistinguishable. 
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